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Abstract 

The high index (441) surface of bismuth has been studied using Scanning Tunnelling Microscopy 
(STM), Angle Resolved Photoemission Spectroscopy (APRES) and spin-resolved ARPES. The 
surface is strongly corrugated, exposing a regular array of (llO)-like terraces. Two surface lo¬ 
calised states are observed, both of which are linearly dispersing in one in-plane direction (kx)^ 
and dispersionless in the orthogonal in-plane direction (/cy), and both of which have a Dirac-like 
crossing at kx=0. Spin ARPES reveals a strong in-plane polarisation, consistent with Rashba-like 
spin-orbit coupling. One state has a strong out-of-plane spin component, which matches with the 
miscut angle, suggesting its possible origin as an edge-state. The electronic structure of Bi(441) 
has significant similarities with topological insulator surface states and is expected to support 
one dimensional Quantum Spin Hall-like coupled spin-charge transport properties with inhibited 
backscattering, without requiring a topological insulator bulk. 
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Bismuth is an interesting material for many reasons [Ml- Being the heaviest stable 
element, spin-orbit interactions are especially significant [S], making it an ideal platform for 
testing fundamental concepts, such as spin non-degeneracy limiting impurity scattering [S] 
and inhibiting a charge density wave [7j. The heavy nuclear weight can also cause a parity 
inversion such that semiconducting bismuth alloys can support an inverted band gap and 
the creation of a topological insulator phase [8] . Such topological insulators have attracted 
much interest [H] and have been suggested as candidates for a range of potential device 
applications dm EH. 

Bismuth and many of its compounds are layered materials. As with graphene, ultra thin 
bismuth is reported to be a two-dimensional (2D) topological insulator [T21 - E5] , supporting 
an edge-localised Quantum Spin Hall (QSH) like state. The vicinal surface of bismuth, 
Bi(ii4) da EH, consists of an array of such edges, and supports a single one-dimensional 
Fermi-contour, which has been revealed to be non-degenerate with respect to spin and is 
reminiscent of a topological insulator state. 

Here, we investigate the surface electronic structure and spin texture of another vicinal 
bismuth surface; Bi(441). A scanning tunnelling microscopy (STM) study was performed, 
and the rhombic surface unit cell of Bi(441) is seen. The (441) surface is observed to 
support a 2 X 2 reconstruction which is well described by a simple truncated bulk, ‘missing- 
row’ model (see Fig. [^a-c)), thus the surface can be viewed as a regular ID array of (110) 
oriented domains, with edges of the (111) bilayers exposed. 

Using angle resolved photoemission spectroscopy (ARPES) we show that the surface 
electronic structure is strongly influenced by the vicinality of the surface. In fact, the Fermi 
surface consists almost entirely of a single ID feature which passes through the surface 
brillouin zone (SBZ) centre, F (see Fig. [^d)). As in the case of Bi(114), the Fermi surface 
shows no sign of the 2x2 surface reconstruction. Whilst the periodicity in is clear, there 
is a lack of any observable dispersion in ky. The lack of periodicity in the surface Fermi 
surface highlights the model ID behaviour of the sample since ky becomes an irrelevant 
parameter in the perfect ID case. 

An additional weak and broad feature (visible as a faint ‘X’) is also present in the ARPES 
measurement (Eig. [^d)). By considering its dispersion with photon energy, it can be at¬ 
tributed to a projected bulk state. This view is strongly supported by tight binding calcula¬ 
tions (see Fig. (e), and the Supplementary Material). There is a lack of bulk symmetry in 
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FIG. 1: (a) Truncated bulk model of the surface, showing the 1x1 (upper half) and the missing 
row (lower half) reconstruction. The 1x1 unit cell is superimposed in green, (b) STM image 
(bias voltage = 0.091 V, tunnelling current = 90 pA.), similarly scaled with the same unit cell 
overlaid, (c) side-view of the 2x2 truncated bulk model and the measured STM height profile, 
(d) Fermi surface map performed with photon energy hu — 70 eV, showing the 1x1 SBZ, and the 
F points. Bright colours (light blue, white, yellow) indicate high intensity and dark colours (black, 
dark red) indicate low intensity, (e) Tight binding calculated constant energy surface, indicating 
where the projected bulk states are close enough to the measured Fermi surface for a weak intensity 
contribution to be expected (using the parameters of Liu and Allen [l8]). 

the ^-direction which hinders a simple identification of = 0. However, the tight binding 
calculated bulk band projection allows F to be identified. 

The single ID feature in the Fermi surface raises the question of the nature of this state. 
Strong spin-orbit coupling is expected [1], and states reminiscent of topological insulator 
surface have been observed on other bismuth [ini uni ED] (and Bia;Sbi_a; [201121]) surfaces. 
However, if bulk bismuth is described as topologically trivial [22H21] then an unpaired Fermi 
level crossing is not expected here. 

Since the surface can be viewed as an array of (111) bilayer edges, an alternative approach 
is to consider the arguments of topology in a similar way; bilayer bismuth (111) is described 
as topologically non-trivial in its 2D bulk, thus its ID edges can be expected to support a 
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topological state. In a simplified view, the high index surfaces of the form 1) behave 
as an array of such (111) bilayer edges, coupled only by van der Waals forces. The lack of 
strong interaction between the layers is manifest in the ID Fermi surface, thus it is perhaps 
reasonable to expect that the ID state inherits the characteristics of the bilayer’s topological 
edge state [HI [25] . 

Measurements made at a range of photon energies (Fig. reveal the existence two 
Dirac-like states, with very similar dispersions in the direction. Their cross-sections vary 
strongly, thus typically only one state dominates the measurement. At particular photon 
energies, both states are visible, for example at hu = 26 eV (Fig. [^c)). A likely cause of this 
strong variation in intensity is a corresponding variation in the matrix element describing 
the photo-excitation to a bulk-like final state |26| . The lack of dispersion with photon energy 
(or equivalently, k± EH) confirms that these states are surface or edge states. Our tight 
binding calculations of the bulk states (which are otherwise in good agreement with the 



- 0.1 0.0 0.1 


k, (A-') 




- 0.1 0.0 0.1 

k, (Ay 



FIG. 2: (a), (b), (c) and (d) ARPES measurements performed at ky^O and photon energies 
18, 22, 26 and 32 eV, respectively, (e) Schematic depicting the presence of ‘inner’ (white) and 
‘outer’ (black) Dirac-like states. 
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FIG. 3: (a) Fermi surface collected with photon energy 19 eV. Yellow horizontal lines indicate 
where the slices shown in (b), (c) and (d) have been extracted (fcy=+0.33, fcy=0 and ky=-0.33 
respectively). 

measurements) do not reproduce these states, further supporting the notion that neither of 
these states are bulk-like. 

The photoemission intensity at the Fermi level (Fig. [^a)), collected at a photon energy 
of 19 eV, shows that this state has a weak dependence on ky. This is best seen in constant ky 
slices taken at ky=+0.33, ky=0 and ky=-0.33 A“^, which reveal a pair of linearly dispersing 
branches (see Fig. [^b-d)), with Dirac-like crossings at A:x=0 with Eb^150, Eb^150 and 
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(a) hv=9 eV 



FIG. 4: (a) ARPES measurement performed at /iz/ = 9 eV showing where two Spin Resolved 
MDGs (SR-MDCl and SR-MDC2) were collected (green dashed lines), (b) and (c) the raw spin 
polarisations collected at the Fermi level, and at = 0.2 eV, respectively, showing the in-plane 
(blue-red) and out-of-plane (grey) spin components. The red axis indicates corrected polarisation^ 
found by subtracting a non-polarised background and correcting for the Sherman function of 17 %. 
(d) ARPES measurement with the reconstructed spin vectors superposed, (e) In-plane (blue- 
red) and out-of-plane (grey) spin components from SR-MDCs atop an ARPES measurement at 
hi/ = 26 eV, where the contribution of the outer state is stronger. In both cases /cy=0. (f) 

Similarly scaled polarisation curves for the four SR-MDCs numbered in (a-c) and (e). (g) Schematic 
representation of the measured spin vectors, and the positions at which they are measured. 
























Eb^O meV for ky=+0.33, ky=0 and ky=-0.33 A“^, respectively. 

In fact, over the range of ky probed (Fig. |^a)), only a very weak dispersion is seen. 
Instead the intensity of the state at deeper binding energy, and the state with the crossing 
nearest the Fermi level, depends on ky. Thus, it appears that both states have a true 
ID nature, and very similar linear dispersions in k^- The band structure of both states is 
presented schematically in Fig. [^e). 

Since strong spin-orbit coupling is expected in this material m El uni [23, we probe the 
spin polarisation of the states observed here by Spin-ARPES, using a Mott-type detector 
[22] • The small k^ separation of the observed states requires good experimental resolutions, 
with measurements at low photon energy being preferable. For example, at photon energy 
9 eV, the spin polarisation can be efficiently measured - the state with its crossing point 
closest to the Fermi level (marked in white in Fig[^e)) is most intense, and hence dominates 
the ARPES measurement, although a very weak remnant of the outermost state is barely 
visible (see Eig. |^a)). Correspondingly, this state dominates the spin-ARPES measurement, 
and hence its polarisation can be found unambiguously. The spin-ARPES measurements 
were performed as two momentum distribution curves (MDGs), one at the Eermi level, and 
the other around 200 meV below. A 2D Mott detector is oriented such that the out-of-plane 
spin component, and the in-plane component in the y-direction are measured. The raw 
spin-polarisation P is given by P = x(Ai+ — N-)/{N+ + N-), where N+ and N- are the 
number of counts recorded by the detectors in the relevant positive and negative directions 
(i.e. ‘up’ and ‘down’ for the in-plane polarisation and ‘left’ and ‘right’ for the out-of-plane 
polarisation). Since the four detectors are not equally sensitive, a sensitivity factor y is 
found by evaluating y = (A^+ -|- A’_)/(Ai+ — N_) well above the Fermi level. In other words, 
the polarisation of the state is evaluated relative to a background measurement. It should 
be noted that no other corrections are made; no Sherman function is assumed, hence ‘raw’ 
polarisation is stated. 

At photon energy 9 eV (Fig. |^a)), two MDGs are selected to cross the state above and 
below its degenerate crossing at k^ = 0. In both spin MDGs, the maximum polarisation is 
seen at k^ ~ ±0.02 (Fig. |^b) and (c)), corresponding to the maximum intensity of the 
state. In both cases, the out-of-plane component (grey) shows no significant polarisation, 
whereas the in-plane component (blue-red) shows a strong sign reversal; negative to positive 
near the Fermi level, and positive to negative at Eb ~ 0.2 eV, typical of a Rashba-like 


coupling. The raw spin polarisation is between 6 and 8 %, which, after allowing for a non¬ 
polarised background, and correcting for the Sherman function (~ 17 %) indicates that the 
state is very strongly polarised. 

In order to observe the spin-polarisation of the second state (with its crossing point at 
around Eb=200 meV and schematically depicted in black in Fig. [^e)), the observed strong 
intensity variation with photon energy is exploited. At hi^ = 26 eV, this state appears more 
intense, and once again spin MDCs can be performed at the Fermi level, and at higher 
binding energy (~ 0.25 eV). The maximum observed polarisation at the Fermi level is now 
seen at ~ 0.05 corresponding to the maximum intensity of the outer state (Fig. 
I^e)). As before, a strong reversal in the in-plane polarisation (from negative to positive at 
the Fermi level) is seen, indicating that the in-plane polarisation of both states is the same. 
Contrary to the previous case, the out-of-plane component is now of a similar magnitude 
to the in-plane component, indicating that the spin vector is approximately 45° to the 
sample surface. The MDC at higher binding energy shows a weaker polarisation reversal, 
probably because of the proximity of a bulk state contributing unpolarised intensity to the 
background. However, this MDC shows a clear reversal of both the in-plane and out-of-plane 
components; from positive (at negative k^) to negative (at positive k^)- 

The results from the spin MDCs are summarised in Fig. |^f). The measured polarisations 
are superposed on the ARPES data for the two relevant photon energies, as well as being 
represented schematically. At the Fermi level, all states have a simple Rashba-like in-plane 
component, whereas the outer branches also have a large out-of-plane component. 

The picture of the electronic structure presented here is broadly consistent with other ID 
systems in which SOC is important |30l |3T] . For example, previous studies of ID metallic 
structures on high index silicon surfaces show Rashba-like coupling and a ID Fermi contour 
ffl. However, these metallic structures on high index silicon surfaces also show a coupling 
of the ID states which is not seen for Bi(441), despite the separation of the ID chains being 
significantly less in the present case |32]. In our case of Bi(441), the Fermi contour shows 
ideal non-coupled ID behaviour. 

In the present case of Bi(441), the observed surface states have a dispersion which is 
more akin to topological insulator surface states than the Rashba coupled free-electron like 
parabola of decorated vicinal silicon [HOI El] • Since alloying bismuth, for example forming 
Bii_a;Sba;, can result in the formation of topological insulator phase [SH|, it is reasonable to 
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assume that Bii_a;Sba:(441) could support an odd number of non-trivial states derived from 
the states observed here. 

The implications of this surface electronic structure for the charge/spin transport prop¬ 
erties are important. The spin texture of the bandstructure is such that a reversal of 
requires an accompanying spin reversal, i.e. e{kx,'t) ^ ^{—kx,-l)- This means that un¬ 
less there is a mechanism for exchanging spin angular momentum with an impurity, the 
probability of a charge carrier (with E ~ Ep and k^ ~ k^^p) being backscattered by an 
impurity approaches zero. Furthermore, there are no other surface states at the Fermi level 
which allow scattering process which conserve |k| and E, thus the probability of momentum 
conserving scattering events is expected to be vanishing small. Indeed, in our STM investi¬ 
gations, there are no visible indications of scattering around impurities. The lack of defect 
scattering, the spin-texture, and the ID Fermi contour suggest that Bi(441) should have 
transport properties akin to quantum spin hall (QSH) materials |31j , which are desirable for 
emergent spintronic applications PES]. 

Although the surface properties of Bi(441) show much promise for spintronics, there is a 
sufficient contribution of unpolarised bulk states to dominate the transport [36] . In order to 
utilise the surface spin-transport properties in a real device, it would be necessary to reduce 
the contribution from the bulk. Two possibilities methods to achieve this are bulk alloying 
(for example with Sb, to open bulk band gap [HZI) and by growth of a vicinal thin layer on 
a semiconducting substrate such as silicon [SHI EH] • 

The existence of a finite out-of-plane spin component is also intriguing. In the simple 
topological insulator case, or indeed for simple Rashba-like SOC, the spin vector at the 
Fermi level is perpendicular to both the surface and the momentum vector, yielding an 
entirely in-plane spin vector. For the high index surface Bi(114), a significant out-of-plane 
component (30°) has been reported [IH], which is comparable to the present case of around 
45° for Bi(441). In both cases, this coincides with the miscut from the (111) plane. In 
other words, although superficially similar to the decorated vicinal silicon case [30] El], the 
two states observed here are not simply SOC branches of the same state, but rather have 
quite separate origins; it is tempting to speculate that the state with a strong out-of-plane 
component originates from a topological-like edge state of a (111) layer [H], and the state 
lacking an out-of-plane component to be the (441) surface state. 

In conclusion, we have presented the electronic structure of a high index surface of 
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bismuth, with a strong ID surface corrugation. Two ID states are observed, both of which 
are non-degenerate with respect to their spin. Bulk bismuth may not be a topological 
insulator, but is nonetheless close to a topological-insulator transition, thus the states 
observed can be rationalised as topological-insulator-like, linearly dispersing, spin-polarised 
states and are expected to support QSH-like coupled spin-charge transport, desirable in 
emergent spintronic applications. 

Acknowledgements: We thank N.A. Vinogradov and A. B. Preobrajenski for facili¬ 
tating the motorisation of the SR-MDCs, and for helpful discussions. S.C. acknowledges 
funding from MAX-IV laboratories and Ph.H. acknowledges the VILLUM fonden for 
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Methods: Samples were prepared by mechanically cutting and polishing a bulk single 
crystal 45° from the (111) natural cleavage plane. Following mechanical polishing, the orien¬ 
tation was confirmed by Laue diffraction. The sample was then electropolished, introduced 
into the ultra-high vacuum chamber and cleaned by multiple cycles of Ar+ sputtering and 
very gentle annealing until contaminants were beneath the detection threshold. 

ARPES measurements were performed at beamline SGM3 of the synchrotron ‘ASTRID’ 
in Aarhus, Denmark m- Low photon energy ARPES and Spin-ARPES measurements 
were performed at beamline 13 of MAX-III [25] • STM measurements were performed at 
Freie Universitat Berlin using a custom-made instrument working in ultrahigh vacuum and 
at low temperature (5 K). Tight Binding calculations were performed using the parameters 
of Liu and Allen [18] . 

Spin-ARPES were measured using linear-vertical polarised light which was incident upon 
the sample at 75° from the surface normal. In other words, the polarisation vector of the 
light is predominantly out-of-plane relative to the sample and with a small component in 
the ky direction. Measurements were also made using the same sample geometry and linear- 
horizontal light (i.e. in the in-plane direction) and the measured electron spin was the 
same within the uncertainty of the experiment. 


* quantum.wells@gmail.com 


11 



[1] V. S. Edel’inan, Physics Letters A 210 , 105 (1996). 

[2] F. Y. Yang, K. Liu, K. Hong, D. H. Reich, P. C. Searson, and C. L. Chien, Science 284 , 1335 
(1999). 

[3] E. 1. Rogacheva, S. N. Grigorov, O. N. Nashchekina, S. Lyubchenko, and M. S. Dresselhaus, 
Applied Physics Letters 82 , 2628 (2003). 

[4] P. Hofmann, Progress in Surface Science 81, 191 (2006). 

[5] X. Gonze, J. P. Michenaud, and J. P. Vigneron, Physica Scripta 37, 785 (1988). 

[6] A. Strozecka, A. Eiguren, and J. 1. Pascual, Phys. Rev. Lett. 107 , 186805 (2011). 

[7] T. K. Kim, J. Wells, C. Kirkegaard, Z. Li, S. V. Hoffmann, J. E. Gayone, 1. Fernandez- 
Torrente, P. Haberle, J. 1. Pascual, K. T. Moore, et ah. Physical Review B 72 , 085440 (2005). 

[8] H. Zhang, C.-X. Liu, X.-L. Qi, X. Dai, Z. Fang, and S.-C. Zhang, Nat Phys 5, 438 (2009). 

[9] J. E. Moore, Nature 464, 194 (2010). 

[10] J. J. Cha, J. R. Williams, D. S. Kong, S. Meister, H. L. Peng, A. J. Bestwick, P. Gallagher, 
D. Goldhaber-Gordon, and Y. Cui, Nano Letters 10, 1076 (2010). 

[11] F. Xiu, L. He, Y. Wang, L. Cheng, L.-T. Chang, M. Lang, G. Huang, X. Kou, Y. Zhou, 
X. Jiang, et ah, Nat Nano 6 , 216 (2011). 

[12] Z. Qiao, W.-K. Tse, H. Jiang, Y. Yao, and Q. Niu, Phys. Rev. Lett. 107, 256801 (2011). 

[13] M. Wada, S. Murakami, F. Freimuth, and G. Bihlmayer, Phys. Rev. B 83 , 121310 (2011). 

[14] S. Murakami, Phys Rev Lett 97 (2006), ISSN 0031-9007. 

[15] L. Chen, Z. F. Wang, and F. Liu, Phys. Rev. B 87, 235420 (2013). 

[16] J. W. Wells, J. H. Dil, F. Meier, J. Lobo-Checa, V. N. Petrov, J. Osterwalder, M. M. Ugeda, 
1. Fernandez-Torrente, J. 1. Pascual, E. D. L. Rienks, et ah, Phys. Rev. Lett. 102 , 096802 
(2009). 

[17] D. Leuenberger, H. Yanagisawa, S. Roth, J. H. Dil, J. W. Wells, P. Hofmann, J. Osterwalder, 
and M. Hengsberger, Phys. Rev. Lett. 110 , 136806 (2013). 

[18] Y. Liu and R. E. Allen, Physical Review B 52, 1566 (1995). 

[19] C. R. Ast and H. Hochst, Physical Review Letters 87 , 177602 (2001). 

[20] S. Agergaard, C. Spndergaard, H. Li, M. B. Nielsen, S. V. Hoffmann, Z. Li, and P. Hofmann, 
New Journal of Physics 3 , 15.1 (2001). 

[21] X.-G. Zhu, M. Stensgaard, L. Barreto, W. S. e Silva, S. Ulstrup, M. Michiardi, M. Bianchi, 
M. Dendzik, and P. Hofmann, New Journal of Physics 15 , 103011 (2013). 


12 



[22] P. Hofmann, J. E. Gayone, G. Bihlmayer, Y. M. Koroteev, and E. V. Chulkov, Physical Review 
B 71 , 195413 (2005). 

[23] J. C. Y. Teo, L. Fu, and C. L. Kane, Physical Review B 78, 045426 (2008). 

[24] Y. Ohtsubo, L. Perfetti, M. O. Goerbig, P. L. Fevre, F. Bertran, and A. Taleb-Ibrahimi, New 
Journal of Physics 15 , 033041 (2013). 

[25] I. K. Drozdov, A. Alexandradinata, S. Jeon, S. Nadj-Perge, H. Ji, R. J. Cava, B. A. Bernevig, 
and A. Yazdani, ArXiv e-prints (2014), 1404.2598. 

[26] J. A. Miwa, P. Hofmann, M. Y. Simmons, and J. W. Wells, Phys. Rev. Lett. 110 , 136801 
(2013). 

[27] F. J. Himpsel, Appl. Opt. 19, 3964 (1980). 

[28] Y. M. Koroteev, G. Bihlmayer, J. E. Gayone, E. V. Chulkov, S. Bliigel, P. M. Echenique, and 
P. Hofmann, Physical Review Letters 93 , 046403 (2004). 

[29] M. H. Berntsen, P. Palmgren, M. Leandersson, A. Hahlin, J. Ahlund, B. Wannberg, M. Mans- 
son, and O. Tjernberg, Review of Scientific Instruments 81, 035104 (2010). 

[30] 1. Barke, F. Zheng, T. K. Riigheimer, and F. J. Himpsel, Phys. Rev. Lett. 97 , 226405 (2006). 

[31] C. Tegenkamp, D. Liikermann, H. Pfniir, B. Slomski, G. Landolt, and J. H. Dil, Phys. Rev. 
Lett. 109 , 266401 (2012). 

[32] 1. Barke, T. K. Rugheimer, F. Zheng, and F. J. Himpsel, Applied Surface Science 254 , 4 
(2007). 

[33] X.-G. Zhu and P. Hofmann, Physical Review B 89, 125402 (2014). 

[34] B. A. Bernevig and S.-C. Zhang, Physical Review Letters 96 , 106802 (2006). 

[35] D. Pesin and A. H. MacDonald, Nat Mater 11, 409 (2012). 

[36] J. W. Wells, K. Handrup, J. F. Kallehauge, P. Bpggild, M. B. Balslev, J. E. Hansen, P. R. E. 
Petersen, and P. Hofmann, Journal of Applied Physics 104 , 053717 (2008). 

[37] F. Nakamura, Y. Kousa, A. A. Taskin, Y. Takeichi, A. Nishide, A. Kakizaki, M. D’Angelo, 
P. Lefevre, F. Bertran, A. Taleb-Ibrahimi, et ah, Phys. Rev. B 84 , 235308 (2011). 

[38] J. W. Wells, J. F. Kallehauge, and P. Hofmann, J. Phys. Cond Matt. 19, 176008 (2007). 

[39] W. Ning, F. Kong, Y. Han, H. Du, J. Yang, M. Tian, and Y. Zhang, Sci. Rep. 4 (2014). 

[40] S. V. Hoffmann, C. Spndergaard, C. Schultz, Z. Li, and P. Hofmann, Nuclear Inst, and Methods 
in Physics Research, A 523, 441 (2004). 

[41] As with the Bi(lll) edge state, the dispersion is approximately linear around the BZ centre. 


13 



and the crossing point at = 0 is close to the Fermi level. However, the band dispersion 
dE/dkx observed here is significantly steeper [13] 


14 


